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Abstract 

We present a simple physical description of the off-diagonal gluon distribution g{x,x',Q'^) 
and splitting function Pgg(x,x') which enter the amplitude for diffractive vector meson 
production. We study the off-diagonal effects both on the evolution in and on the in- 
put distribution. We predict the ratio R of the off-diagonal to the diagonal gluon density, 
x'g{x,x',Q'^)/xg{x,Q'^), as a function of the kinematic variables. 
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1. Introduction 

It is well known that the cross section of hard scattering processes (such as deep inelastic 
scattering, the production of large px jets, etc.) can be written as the sum of parton distribu- 
tions multiplied by the cross sections of hard subprocesses calculated at the parton level using 
perturbative QCD. That is we can factor off the long distance (non-perturbativc) effects into 
universal, process independent, parton distributions specific to the incoming hadron. 

Strictly speaking a parton distribution is the diagonal element of a density matrix, or the 
diagonal element of an operator O in the Wilson operator product expansion (OPE). To be 
more specific let us start with the Fock decomposition of the wave function of an incoming 
hadron, say a proton 



where k — qqq, qqqg, . . . and n is the number of partons in the Fock state ijjk- Each parton is 
described by a vector with coordinates which specify its momentum, colour, spin, flavour 
etc. To calculate an inclusive cross section we must sum over the coordinates of all the partons 
except for the one which participates in the hard subprocess. Thus, for example, the gluon 
distribution is the diagonal element 



where as usual the operator product GG denotes the product of an annihilation operator G and a 
creation operator G^. The product "counts" the number of gluons in the Fock states. Contrary 
to the conventional notation the operator G corresponds to a gluon with fixed momentum 
fraction x, not to a point r in coordinate space. 

On the other hand there are a set of reactions which are described by the off-diagonal 
elements of the density matrix, say {N'\0\N) where the momentum, helicity and/or charge of 
the outgoing nucleon N' are not the same as those of the ingoing nucleon N. Examples are : 

(a) virtual photon Compton scattering (7*^ 7p) or vector particle electroproduction 
(7*p — i> Vp) — examples of the latter processes include ■j*p — > pp, ■y*p — > Zp, 7*p —>■ W^n 
and 7*p — >• {ii^ii~)p; 

(b) processes which involve the BFKL amplitude at large \t\ — \p — p'P 7^ 0; 

(c) the polarized structure function g2, or alternatively the function g± = gi + g2 for trans- 
versely polarised nucleons. 

In each case we have to study the evolution of the matrix element of the two gluon operator 
GG (or the two quark operator qq) evaluated between non-diagonal states, ^ a[. In case 
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(a), 7*p — s> Vp, the important difference in {p'\GG\p) is between the longitudinal components 
of the incoming and outgoing proton momentum, where 

tmin = {P-P)\\ = - ( — 1 m^. (3) 

is the virtuality of the photon, W is the 7*p centre of mass energy and My is the mass of 
the produced vector meson. On the other hand in case (b) the large difference is due to the 
transverse components 

|t|^(p-p')T, (4) 
whereas in case (c) the difference is due to the nucleon helicities, (A'|gg|A) with A 7^ A'. 

To describe these processes we need to consider off-diagonal elements of the density matrix 

„ n 

{N'\d\N) = E y ■■■,<) OMai, • • • , «n) n <^(«* - <) d< (5) 

k ij^g 

where the coordinates of the Fock state and ip* are all equal (a^ = a'j) except for the parton, 
say i = g, which participates in the hard scattering. A topical example is the off-diagonal gluon 
distribution g{x,x') = {p'\GG\p) for which the longitudinal momentum fraction carried by the 
gluon {i = g) differs on the \eit-{ip*) and right-('?/') hand sides of the amplitude, Xg 7^ x'g. 

In this paper we concentrate on just such off-diagonal parton effects which occur in the 
description of diffractive vector meson production, since these processes appear to be the most 
accessible experimentally. We introduce the appropriate variables in Section 2. Off-diagonal 
effects can change both the evolution and the input distributions. We study off-diagonal evolu- 
tion in Section 3. Although the off-diagonal splitting functions are already known, this section 
contains a particularly simple and physical derivation of Pgg{x,x'). The off-diagonal effects 
on the starting distributions (which are used for the evolution) are discussed in Section 4. In 
Section 5 we present numerical estimates of the size of these off-diagonal effects as functions of 
X, x' and . Finally in Section 6 we present our main conclusions. 

2. Diffractive electroproduction of vector mesons 

Perhaps the most relevant reactions, at present, are the quasielastic diffractive electropro- 
duction of vector mesons, 7*p — > Vp. These processes are now being studied with increasing 
precision at HERA. It has been argued 0] that the amplitudes for these reactions are, to a 
good approximation, proportional to the gluon distribution, g{x,Q ), that is the conventional 
(diagonal) parton distribution measured in the global analyses of data for deep inelastic and 
related hard scattering processes. To lowest order the 7*p Vp amplitude can be factored into 
the product of the 7* — > qq transition, the diffractive scattering of the qq system on the proton 
via two-gluon exchange, and finally the formation of the vector meson V from the outgoing 
qq pair. The crucial observation is that, at high cm. energy W , the scattering on the proton 
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occurs over a much shorter timescale than the 'J ^ qq fluctuation or V meson formation times. 
Observation of YP ~^ at cm. energy W therefore probes the gluon distribution g{x, Q ) 
with 



g2 + 



X = Q' = z{l-z)Q'' + h;'r + ml (6) 

where z, and 1 — z, —kx are the g, q longitudinal momenta fractions and transverse momenta 
with respect to the incoming photon, and is the mass of the quark. Indeed if this simple 
model can be theoretically substantiated, the dependence of the cross section on the square of 
the gluon density means that data for 7*p Vp (and for the photoproduction of heavy vector 
mesons e.g. 7p J/'ipp) offer an especially sensitive probe of the gluon. 

However we have noted that 7*p — > Vp is actually described by a non-diagonal distribution 
g{x,x') and so it is important to quantify^ the difference between this improved description 
and the approximation of using the conventional, diagonal, gluon distribution. To be specific 
let us study 7*p —>■ pp at high Q^. Due to the proton form factor, the transverse momentum of 
the p meson is strongly limited (p^ <^ Q^) and the difference between the amplitude and the 
one with diagonal gluons arises predominantly from the transfer of the longitudinal momentum 
which is needed to put the p meson on-mass-shell. The partonic decomposition of the process 
is sketched in Fig. 1. Also shown are the relevant kinematic variables. 

Consider the imaginary part of the amplitude, which is the dominant component. Then the 
intermediate particles are on-mass-shell as indicated by the crosses on the lines in Fig. 1. In 
order to put the quarks {k and k) on-shell we require the fraction of the proton's longitudinal 
momentum that is carried by the gluon to be 

^ = — TJT^- (7) 



Mqq, the mass of the qq pair, is given by 



Mj, = (k + ky ^ -^—^ + ^--^ - 4 (8) 

where recall that {z, kx) and (1 — z,kT = It — kx) specify the g, q kinematics with respect to 
the photon. 

On the other hand for the qq —>■ p transition the difference between the mass of the vector 
meson My and the qq mass M^g is very small. Now in the leading log approximation (LLA) 
the transverse momenta are strongly ordered, and in particular <^ Q^. Therefore 

x' = « X. (9) 



Thus we have a longitudinal momentum difference 

Q: 



5x = Xi-Xi ^ — — — , (10) 



^Recent work on this subject can also be found in ref. 
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which persists along the whole ladder, see Fig. 1. 

The variable x' is not directly measurable. Since the mass of the qq system, M^g, is not 
fixed, the value of x' is smeared out by the integration over the quark momenta. It means that 
observation of diffractive vector meson production, 7*]? Vp, measures the off-diagonal gluon 
distribution over some interval of x'. However this is not such a disadvantage as it seems. To 
leading log accuracy the range of x' is much less than x and we may safely put x' = (see 0). 

We wish to study the InQ^ evolution of the off-diagonal parton distributions indicated by 
the ladder structure in Fig. 1. First we note that this is part of a more general evolution. In fact 
for the electroproduction of a massive vector boson, say 7*p —>■ Zp, the off-diagonal evolution 
equation [|, |^, ||, |^, || describes not only the deep inelastic evolution for spacelike Q^, which is 
our concern here, but also the Brodsky-Lepage leading In My evolution of the vector boson 
wave function which occurs in the timelike region.| For instance for heavy vector boson 
production 7*p — > Vp we have evolution arising from the sum of (a^ In My)" terms in the 
region of x' < 0. Since now both the x and x' gluons go in the same direction, we are in the 
timelike region. However, we lose at least one InQ^ at the cell where x' changes sign. When x' 
goes negative the parton reverses direction. Therefore there are no coUinear logs at this point 
and, moreover, x' < kinematics is inconsistent with causality in the LLA. Indeed the x' < 
domain corresponds to the decay of the heavy vector boson into the virtual photon and two 
gluons, for example Z 1*99- In the LLA it is described by the point-like production of two 
massless (on-shell) gluons. However these two gluons cannot interact through s channel gluon 
exchange (as shown in Fig. 1): the x' < gluon has insufficient time to overtake the gluon 
X. A more formal and detailed discussion is given by Radyushkin [^]. So at leading order we 
must terminate the evolution at x' ~ and the amplitude of the process sketched in Fig. 1 is 
proportional to the off-diagonal gluon distribution^ x'9{x,x') with 

^ - — — ' ^ - ° 

In the evolution we keep the full x dependence. For the process under consideration, 
1*P ~^ PP ci't HERA energies, the values of x are very small, x ^ 10~^. In this small x region the 
leading contribution has a double logarithmic form in which at each iteration we have not only 

^This is analogous to the Gribov-Lipatov relation which connects the splitting functions for ep eX 
deep inelastic evolution in the spacelike region with the splitting functions which occur in e"'"e~ annihilation 
(or, to be more precise, in e^e^ qX) at moderate x in the timelike region. 

^We use a factor x', rather than x, to avoid the 1/x' singularity for x' ^ x. It is informative to amplify 
this comment. Recall that for the diagonal gluon density we use xg{x) to avoid the Bremsstrahlung singularity, 
where g{x)dx describes the probability to find a gluon in an interval dx of momentum space. When we express 
observables in terms of the gluon distribution we find that the relevant region of the dx integration is of the 
order of a;. It would therefore have been more natural to introduce parton densities xg{x) in the place of g{x). In 
the off-diagonal case, the region of the dx integration (which corresponds to dag in (|^)) is limited by minimum 
{x,x'} = x' . So the appropriate density is x'g{x,x'). 
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In but also a In 1/x arising from the integration over longitudinal momenta. This double 
leading logarithmic part of the amplitude,^ 

x'g{x,x') ~ Cn(«5ln- InQ^)", (12) 

corresponds to the domain in which the longitudinal, as well as transverse, components are 
strongly ordered, Xi -C Xi_\. As in this domain Xj ~ x,, there is no x' dependence and the 
off-diagonal distribution becomes the diagonal one.0 

Therefore we can consider the use of the diagonal parton density as a zero approximation 
and study the difference between the diagonal [bx — > 0) and off-diagonal amplitudes. That is 
we will calculate the ratio 

R = — ' ^ where x' = x — Sx. (13) 

xg{x) 

There are two possible sources of difference between the diagonal and off-diagonal densities 
which causes the ratio R to differ from unity. First the input distributions at the starting 
scale Qq of the evolution may differ and, second, there is a difference in the evolution itself on 
account of the difference in the splitting functions. It is convenient to discuss the effect on the 
evolution first. 



3. OfF-diagonal evolution equation and splitting function 



The off-diagonal splitting functions which describe the log evolution of partons have been 
widely discussed for example [|12|, ^ H, |, |^. They have been calculated by the QCD-string 
operator approach |^, ^ or by using the conventional LLA in which the partons are put on- 
shell, that is by the QCD extension of the Weizsacker- Williams technique^. So the off-diagonal 
splitting functions P{z, z') are known. 

Since the gluons play a dominant role in the small x HERA domain it should be a reasonable 
simplification to consider pure gluon evolution. The evolution equation for the off-diagonal 
gluon density is then 



dg{x,x') OisiQ"^) dz 



[ -P,,{z,z')g{x/z,x'/z') (14) 
Jo z 



^Recall that the double logs come from the integration over the rapidity of the s channel gluon (and from 
the integral over the strongly ordered transverse momenta). In the notation of ref. |^ the upper limit in the 
dxs/xs integration is given by the largest of the t channel momentum fractions in the i"' cell, that is by Xi 
rather than by x^. Hence we find In l/x and not In 1/x' . 

^This was the conclusion of the first paper on off-diagonal distributions. 

^The splitting functions were originally discussed in ||l3|l , but unfortunately there were mistakes in the 
formulae. 
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where in the notation of Fig. 1 



z' = 4^ = ^ . (15) 



Note that in (0) there is no integration over z' . To see why this is so recall that the dz' j z' 
integration is equivalent to the dx' j x' integration and that x' = x — bx. Since bx is fixed by 
([l0|) then z' is fixed also. The splitting kernel is given by [T^, |^ 



/N / ,N ( 1 1 Z — z' z' — Z 

Pgg{z,z') = 2iVcz(l-/) 1 + — + ———-- + — —— + 



{l-z){l-z') 2z'{l-z) 2z{l-z')} 

- {V + V')6{1- z). (16) 

The integrals V and V in the virtual term correspond to the left- and right-hand sides of the 
amplitude. We have 



V = Nc I zdz 







^ 1 

z(l-z) + - + 



■ , (17) 

z 1 — z] 

In an attempt to give physical insight of the structure of Pgg we will, below, present a particu- 
larly simple derivation of (p^G]). 

First we notice that if z = z' then Pgg{z,z') reduces to the conventional splitting function 
Pgg{z). Recall that in the conventional (diagonal) case, the 1/(1 — z) singularity as 2 — 1 in 
the real part of the kernel (which corresponds to the emission of a soft gluon with momentum 
fraction tji = Xi{l ~ z)) is cancelled by the singularity in the virtual contribution which is 
needed to restore the normalization of the parton wave function. To be explicit, the real- virtual 
cancellation is _ 

(18) 

which gives a combination that is regular as z — 1. For the gluon C2{G) = Nc = 3. In 
accordance with the conventional Bloch-Nordsieck procedure we must cut off the contribution 
due to soft gluons with small energies, say with E < Eq, using the same cut-off Eo in the real 
and virtual parts of the splitting functions. 

At first sight it appears from (|l^) that for off-diagonal evolution only the pole 1/(1 — z) 
occurs, but not a 1/(1 — z') singularity. On the other hand we see that the normalisation is 
2Nc, which reflects the fact that for fixed 6x = Xi — x'^ and z ^ 1 (that ) the 

value of z' = x'Jx[_i also goes to 1. That is the 1/(1 — z) pole represents both the 2; — > 1 and 
2;' — > 1 singularities. So the total soft gluon singularity in the real part with a factor 2Nc is 
balanced by the sum of the V and V' virtual terms. However care is needed. We note that 
in the off-diagonal case we have different momentum fractions in the "left" and "right" virtual 
terms, that is x x'. Thus we must use slightly different cuts for 

Zi = Xi/xi-i and z[ = x'Jx[_i. (19) 
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Suppose that we use 1 — Zi > A in the virtual term V in the left amplitude, where A corresponds 
to the same bound E > Eq that is imposed on the real soft gluon emissions, then we must use 
a cut 

in the virtual term V in the amplitude on the right. 

It is instructive to sketch the derivation of the off-diagonal splitting function Pgg{z,z') of 



(16). The relevant triple gluon vertex is shown in Fig. 2. To leading order accuracy we can put 



the gluons on-mass-shell and then the vertex function 

^pup = 8^v{p + k)p + 5up{p-2k)^ + 5pp{k-2p)^ (21) 

reduces to simple forms depending on the directions of the polarisation vectors. For each gluon 
we take one of the two independent polarisation vectors to be e*^"^ in the direction of p x fc. 
Then the other polarisation vector is taken, for the s channel gluon to be e^*) in the direction 
e*^") X q, and for the t channel gluons to be e^*^ and e^*) in the directions e*^") x p and e*^") x k 
respectively. It is convenient to introduce 



abc 



r..p eif) e(^) e^;^ (22) 



where a, b, c must be either n or t to indicate the polarisations of the three gluons. Then the 
only non-zero elements of Tabc are 

Tnnt = Tut = 2kT{-^ + --l) (23) 

1 — z \1 — z z / 

where z = k\\/p\\ = Xi/xi^i. We may make two observations. First we note that ( pB]) satisfies 
the leading logarithm approximation requirement that the helicity factor |rp be proportional 
to k^. Second, it is straightforward to see from (^) that a non-zero Tabc must have an even 
number of n subscripts. We may use (^3]) to obtain the real gluon emission contribution to Pgg. 
We sum over the final, and average over the initial, polarisations, and find 

i yr(z')T(z) = 44/1 + — + - — ^ ^ + ^~ ^' , + f ~ \ I . (24) 

2^ v;v; T| ^ zz' {l-z){l-z') 2z'{l-z) 2z{l-z')j ^ ' 

Finally to obtain the contribution to Pgg given in (|TB]) we must include the kinematical factor, 
\z{l — z'), which is appropriate when writing the evolution equation ([l^) as an integral over 
dz/z. In the limit z' ^ z we see that Pgg{z,z') reduces to the familiar (diagonal) splitting 
function Pgg[z), and that the last two terms in ( p4|) vanish. These terms come from the cross 
products in ri^^{z') Tttt{z). 
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4. Off-diagonal effects on the input distributions 



Suppose that in the diagonal approximation we provide the input at some (low) scale Qq. 
Now the off-diagonal matrix element can be expressed in terms of the Mandelstam variables 
(s, u, t) and the masses or virtualities of the incoming and outgoing particles. It corresponds 
to the elastic parton-nucleon amplitude, say A{s, t, Qq, Qq ), which describes the subprocess 
shown at the bottom of the ladder diagram of Fig. 1. We have to study two effects. One due 
to t 7^ and the second due to the difference in the virtualities of the two outgoing gluons 

SQ' = Q'o' -Ql^Ql f -Ql = -Ql ^, (25) 

\ Xq J Xq 

in the notation of Fig. 1. We consider x <^1. Of course there are important off-diagonal effects 
at large t which originate first from the form factor of the target proton and second from the 
BFKL ln(l/x) evolution at non-zero p^. 

Here we are not concerned with large p'rp vector meson production. However even for zero- 
angle scattering we still have 

\t\ = |tmin| = ml{5xf ^ 0, (26) 

so we have to consider a potential i 7^ off-diagonal effect. Now the scale which controls the 
t behaviour of the amplitude is given by the nearest singularity {t — 4m^) or by the slope 
S in t of elastic parton-proton differential cross section da/dt oc cxp{Bt), where the value of 
B is less than the slope of the elastic pp cross section B < Bpp ~ lOGeV^^. Thus for HERA 
energies, where Sx < 10~^ — 10~^, the effect coming from |tmin| ^ 10~^GeV^ is negligible. 

On the other hand there may be some suppression of the input distribution due to SQ^ ^ 0. 
So far the current "off-diagonal" models of the nucleon would imply a weak dependence on 
Sx = X — x', for 6x <^ 1, at the initial scale, that is before the evolution is taken into account. 
Thus the calculations have taken the input ratio 

x'g(x,x' ,Ql) , . 

at t — train- Howcvcr it is not clear at which value of Ql we should take it!o = 1- The evolution 
causes the ratio R to increase, particularly for x ~ 5x (that is ~ 0). Thus, strictly speaking, 
after we have chosen the starting scale Q^, we should determine Ro{x, 5x) by fitting to the data. 
In other words off-diagonal evolution requires a two- variable input function x/g{x, x', Qq) at an 
initial scale Qq, which originates from the non-perturbative domain. Clearly Ro{x, Sx) depends 
on our choice of Qq. 

Lacking the relevant data to determine Rq{x,Sx), the best that we can do at present is 
to start the evolution from rather small Qq where nevertheless we may use perturbative QCD 
evolution. At these values of Q^ we may expect i?o — 1- The reason is as follows. On the one 
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hand it appears unreasonable to have i? < 1 at small Qq, as in conventional phenomenology the 
elastic parton (gluon) - proton amplitude A{s, t, Q^, Q'q) generally increases as the virtuality 
Qq decreases. On the other hand we may use the Schwarz inequality 

j ij*{x',ai) ip{x,ai) dSg Y[dai < ^ J 0^)^ + \t/j{x, a,)\^^ dag Y[ dui (28) 

to put an upper limit i?o(max) on Rq. We are using the notation of (j^), but for simplicity we do 
not show the sum k over the Fock states. By dcLg we mean integration over all the coordinates 
of the gluon except for dxdx' . The Schwarz inequality (^) gives the upper bound 

^ - - 2xgMl) ■ ^ ^ 

We notice that in (^), written in terms of Fock wave functions, the dependence is not 
explicit. The conventional gluon distribution g{x^Q'^) represents the number of gluons with 
transverse momentum satisfying k"^ < Q^. Thus the argument of plays the role of upper 
limit of the 11 dk'^^ integration in the subspace of the integration over 11 dai. That is the 
evolution scale is controlled by the largest transverse momentum of the partons and not by 
the virtuality Q'^. To very good accuracy we find Roi'oiax.) = 1, throughout our x interval of 
interest, using GRV partons [|l|] at Ql = 0.4GeV^ whereas for MRS(R2) partons H we find 
i?o(niax) ~ 1 in the range = 1.3 — 1.4GeV^ as 5x varies from 10~^ to 10~^. 

Finally, for completeness, we note that the input for the off-diagonal valence quark distri- 
bution has recently been discussed in ref. P|. 



5. OfF-diagonal effects from evolution 

To indicate the effect of using off-diagonal distributions we calculate the ratio 

xg{x,x,Q ) i"3n\ 
^ - xgix,Q^) 

for fixed values of 6x = x — x', using both the leading order off-diagonal and the conventional 
(diagonal) evolution equations for the gluon. We use the gluon of GRV JTSf as input at Ql = 



l.SGeV^. The results for R{x,Sx,Q'^) are shown in Fig. 3 for four values of Sx, in each case 
showing _R as a function of x for = 4, 20 and lOOGeV^. 

The main features of the results shown in Fig. 3 are as follows. First, the ratio R only 
noticeably differs from 1 at small x, close to 6x. For x ^ 10 6x we see that R is always less 
than 1.05. Second, as may be expected, R increases as we evolve up in Q^. To see the reason 
for the increase we must consider the real and virtual contributions to Pgg. The diagonal and 
off-diagonal splitting functions have the same 1/z behaviour and thus lead to the same double 
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log behaviour, a^liiQ^ Inl/x, of the real gluon emission contribution. On the other hand, if 
xg ~ a;~^ {Q'^V ) then the anomalous dimension 

7(A) = g P,,{z) dz (31) 

due to the 'diagonal' Pgg is less than that due to the 'off-diagonal' Pgg, mainly due to the larger 
cut-off A' in the virtual V term as compared to A, see (pOl) . We therefore have a smaller 
remainder in the negative V contribution, after the cancellation of the singularities as in ([ISD . 
This in turn leads to a positive difference in 7(A) and to a faster off-diagonal evolution. 

The third effect apparent in Fig. 3 is the increase of R with 5x. Why is this? The difference 
between the splitting functions AP = Pgg{z, z')—Pgg{z) decreases with decreasing z—z' ^ 6x/x, 
see (p!5| ) and (p!6|). Thus most of the difference comes from the region of evolution where 6x ~ x. 
If we evolve from a value of x = xq which is much larger than 6x, then it is only the last few 
evolution iterations which have 5x ~ x and so contribute io R ^ 1. From this point of view 
we expect that R depends on 5x/x and not on the absolute value of 5x. This is true for very 
small 5x, and hints of this behaviour can be seen by comparing the third and fourth plots of 
Fig. 3, that is those for 5x = 10~^ and 10~^ respectively. However for larger values of 6x, say 
6x ~ 10~^ the evolution starts from xq {xq ~ 0.1) which is already sufficiently close to 6x so 
that all intervals of the evolution contribute to i? 7^ 1, giving a larger value of R. 

Fig. 4 shows the effect of taking the input Ro = 1 at a considerably higher starting scale, 
namely Ql = 4GeV^. The lower two curves in each plot compare the values of R{x, 5x, Q"^) 
obtained using this Ql (dot-dashed curves) with our preferred predictions based on Qq = 
l.SGeV^ input (continuous curves). From the previous discussion we expect the Ql = l.SGeV^ 
input to yield a larger value of R simply because of the longer interval of evolution. However 
we see that the difference disappears during evolution mainly due to the fact that for very large 

the essential contribution to the gluon density comes from the region of initial Xq ^ 6x 
where the off-diagonal effects in the initial conditions are negligible. For interest, we also 
started the evolution from the lowest possible value of Ql, namely Ql = 0.4GeV^, again using 
the GRV gluon. The predictions are shown by the dashed curves, which by chance coincide, 
for an appreciable range of x, with the upper limits which we discuss below. 

The remaining curve in Fig. 4 is obtained from the upper limit (pOj) given by the Schwarz 
inequality. Of course we may use the inequality at any in the perturbative QCD region. 
Clearly calculating the bound at the lowest reasonable will give the tightest constraint on 
R. For this reason we use GRV gluons to calculate i?(max) of ( pOf ) at = l.SGeV^ and then 
evolve both g{x,x') and g{x) to give the "evolved" upper limit shown at the higher values 
in Fig. 4. This procedure is valid because the evolution equations are linear in g. For small 

the difference between R obtained by evolving from the Rq = 1 at Ql = l.SGeV^ (and also 
Ql = 4GeV2) and the upper limit is appreciable. However as we evolve to large the system 
depends less and less on the initial conditions and the predictions approach each other. 

In Fig. 5 we explore the sensitivity of the predictions for R to the choice of the input gluon 
distribution. We compare the predictions based on using the gluon from the MRS(R2) set of 
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partons [16] with those previously obtained using the GRV gluon |15]. In each case we start the 



evolution from Qq = l.SGeV^. The predictions are shown by continuous curves, together with 
the upper limits obtained by evolving from the Schwarz inequality evaluated at Qq = l.SGeV^. 

We may make two observations concerning Fig. 5. First, we see that the MRS and GRV 
predictions for the ratio essentially coincide at = 4GeV^. However the individual off-diagonal 
and diagonal gluons are larger and steeper for GRV than MRS. Since xg ~ has a large 
effective A for GRV than MRS it gives, on average, a larger 6x/x and hence a more rapid off- 
diagonal evolution. Thus the GRV prediction for R lies above that of MRS at = lOOGeV^. 
Second, we see that the MRS predictions lie just below the upper limit, whereas the upper limit 
for GRV is considerably higher. Again this is directly attributable to the steepness in x of the 
GRV gluons. In fact the GRV gluons |]T3[ are "steeper" than that required to fit the new high 



precision measurements of F2 at small x. Therefore the upper limit based on the GRV gluon 
should be regarded as too high. 

Qualitatively our results for R have similar features to those given in ||^. However our pre- 
dictions are appreciably smaller. For example, at x = 1.1 x 10~^, x' = 10^^ and ~ lOOGeV^ 
we obtain R = 1.28 as compared to R = 1.47 of ref. [^. Our study of the off-diagonal effects 
of the input indicates that this is not the origin of the difference.]] 

6. Conclusions 



We have studied the effect of using an off-diagonal gluon distribution g{x, x', Q^) to de- 
scribe diffractive vector meson production. As compared to the conventional diagonal gluon 
distribution we find 

depending on the kinematic variables, see Fig. 3. In particular, R increases with and 
with 6x = X — x' . The enhancement is important phenomenologically since the cross section 
for diffractive vector meson production depends on the square of the gluon distribution. For 
example for J/ip photoproduction, where the effective = Mj^^/A ~ 2.5GeV^ 0, there is an 
enhancement by a factor 1.2 using off-diagonal evolution from = l.SGeV^ with either the 
GRV |15| or the MRS(R2) [T^ gluon as input. This is mainly a normalization effect and does 
not change the dependence of the cross section on the jp centre-of-mass energy W (or on x), 
since R mainly depends on the ratio 6x/x. 

In the absence of data, in Section 5, we studied a range of possibilities for the input distribu- 
tion g{x, x', Qq). We find that the choice of input has little effect on the value of R at large Q'^, 
but that it is not negligible for ~ lOGeV^. We therefore await data for the determination of 
the input distribution. From this point of view the measurement of J/ ip production in a wider 

^We have an even larger difference, both numerically and in the x behaviour, with the results shown in Fig. 2 
of ref. il. 
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kinematic range is particularly attractive. As there appears to be little intrinsic charm in the 
proton, we may hope to use these data at a wide range of W to fix the gluon input over an 
appreciable range of x. 
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Fig.l 

Figure 1: A diagrammatic representation of diffractive p meson production, 7*p — > pp, via a 
two-gluon exchange ladder. The crosses indicate that the particles are on-mass-shell in the 
calculation of the imaginary part of the amplitude. 
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Fig. 2 

Figure 2: The variables of the three-gluon vertex used in the vertex function of (|21|). 
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Figure 3: The ratio R of the off-diagonal to diagonal gluon density as a function of x for four 
values of 6x = X — x' and three values of Q^, obtained from evolving from the GRV gluon ^5 
with RQ = latQl = l.SGeVl 
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Figure 4: The dependence of the ratio of the off-diagonal to diagonal gluon density on the 
input distribution. The dashed, continuous and dot-dashed curves are obtained by evolving 
from Rq = 1 at Ql = 0.4, 1.5 and 4GeV^, respectively, using the GRV gluon |T^. The 
upper limits (represented by dots) are obtained by evolving from Ro(ma.x) of (|29| ) evaluated at 
Qo = l.SGeV^. The predictions corresponding to Qq = 4GeV^ are shown only for comparison 
and should not be included in a realistic estimate of the uncertainty in R. 
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Figure 5: The predictions, together with the upper hmits, based on evolution from the MRS(R2) 
gluon [T^ compared with those obtained by evolving from the GRV gluon [T^. In each case 
we choose Ql = l.SGeV^. 
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